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Abstract 

Wc  consider  Laplace's  equation  in  a  polygonal  domain  together  with  the  boundary 
conditions  that  along  each  side,  the  derivative  in  the  direction  at  a  specified  oblique  angle 
from  the  normal  should  be  zero.  First  we  prove  that  solutions  to  this  problem  can  always 
be  constructed  by  taking  the  real  part  of  an  analjrtic  function  that  maps  the  domain  onto 
another  region  with  straight  sides  oriented  according  to  the  angles  given  in  the  boundary 
conditions.  Then  we  show  that  this  procedure  can  be  carried  out  successfully  in  practice 
by  the  numerical  calculation  of  Schwarz-Christoffel  transformations.  The  method  is  illus- 
trated by  application  to  a  Hall  effect  problem  in  electronics,  and  to  a  reflected  Brownian 
motion  problem  motivated  by  qucucing  theory. 
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1.  The  oblique  derivative  problem  and  conformal  mapping 

Let  n  be  a  polygonal  domain  in  the  complex  plane  C,  by  which  we  mean  a  possibly 
unbounded  simply  connected  open  subset  of  C  whose  boundary  dCi  consists  of  a  finite 
number  of  straight  lines,  rays,  and  Une  segments.  Let  these  sides  be  denoted  by 
Tj,  .  .  .  ,r„,  n>l,  where  T^  is  the  complex  open  interval  (24,2;+,),  for  some  sequence  of 
vertices  z^  €  CU{a=},  and  let  (I  he  everywhere  to  the  left  of  dil  as  dil  is  traversed  in  the 
order  Tj.r^,  •  •  •  .  For  convenience  we  set  zq^z„  and  z„+^^z^.  We  permit  geometries  in 
which  one  or  more  slits  are  cut  into  fi  by  viewing  the  two  sides  of  such  a  sUt  as  distinct 
boundary  segments,  e.g.  T^  and  r^+i.  The  closure  O  is  then  a  manifold  that  cannot  be 
embedded  in  the  plane  (although  it  can  be  reduced  to  an  embeddable  region  by  a  confor- 
mal map),  and  a  function  in  C(fi)  will  have  distinct  values  on  T^  and  r^+i  corresponding 
to  Umits  on  the  slit  taken  from  different  sides.  One  could  also  go  further  and  let  fl  be  a 
polygonal  Riemann  surface,  but  we  will  not  do  this. 

Let  a  set  of  real  numbers  Sj e„  be  given.   At  any  point  z  e  r^,  let  u„{z)  =  — (z) 

dn 

denote  the  inward  normal  derivative  of  a  function  u,  and  u^  =  —  its  tangential  deriva- 
tive along  r^  in  the  positive  (counterclockwise)  sense  with  respect  to  ft.  This  paper  is  con- 
cerned with  the  following  oblique  derivative  problem: 

Problem  O.   Find  functions  u  d  C^(ft)nC^(ftU  Ui^i)  satisfying  Laplace's  equation 

A«(z)  =  0,        z€ft  (1.1) 

together  with  the  homogeneous  oblique  derivative  boundary  conditions 

cos9iM„(z)  +  sin9iU„(z)  =  0,        r  e  T^.  (1.2) 

The  situation  is  illustrated  in  Figure  la.  If  64  =  0  (modir),  we  have  a  Neumann  condition 
on  side  r^,  while  6^  ^  ^(mod'iT)  gives  a  tangential  condition,  which  we  do  not  rule  out. 

Obviously  any  constant  function  is  a  solution  to  Problem  O.  K  ft  were  a  smooth 
domain  with  a  continuous  single-valued  obliquity  function  6(ct)  for  a  e  aft,  the  constants 
would  be  the  only  solutions,  as  can  be  shown  by  consideration  of  Ricmann-Hilbcrt  prob- 
lems [12,13,20].  More  generally,  if  ft  were  smooth  and  e(a)  changed  continuously  by  ATtt 
as  a  traversed  aft,  there  would  be  a  solution  space  of  finite  dimension  raax{l,A'}.  But 
Problem  O  is  quite  different  from  the  analogous  problem  on  a  smooth  domain,  for  no 
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Figure  1.   Problem  O  and  its  solution  by  conformal  mapping.   The  derivative  of  u  at 
the  boundary  of  il  in  the  direction  given  by  the  arrow  must  be  zero. 

boundary  conditions  have  been  specified  at  the  vertices,  nor  have  we  required  regularity 
or  boundcdness  there.  As  a  result  the  space  of  solutions  has  infinite  dimension,  so  that  in 
any  particular  application,  additional  conditions  will  be  needed  to  ensure  uniqueness.  We 
have  intentionally  omitted  such  conditions  from  the  statement  of  Problem  O  for  the  rea- 
son that  their  most  natural  form  varies  from  problem  to  problem.  The  statements  of 
Problems  O^  and  Og  in  Sections  2  and  3  below  incorporate  appropriate  additional  condi- 
tions for  the  two  applications  considered  there. 

The  purpose  of  this  paper  is  to  construct  nonconstant  solutions  to  Problem  O  by  the 
following  process  of  conformal  mapping,  illustrated  in  Figure  1.  Suppose  an  analytic 
function  /  in  HUUil^i  is  found  that  maps  CL  onto  another  polygonal  domain  f(Cl).  (/ 
may  fail  to  be  globally  or  locally  one-to-one,  in  which  case  we  view  f(fl)  not  as  a  subset  of 
C  but  as  a  Ricmann  surface.)  Suppose  further  that  each  image  /(P^)  consists  of  one  or 
more  sides  of  /(H)  oriented  at  an  angle  9^  (modir)  counterclockwise  from  the  real  axis. 
For  each  point  z  e  T^,  either  /'(z)  =  0,  or  /  is  conformal  at  z.  Assuming  the  latter,  we  sec 
that  /  maps  any  curve  that  meets  T^  at  z  at  an  angle  6^  clockwise  from  the  inward  normal 
(Figure  la)  onto  a  curve  that  meets  /(PJ  at  /(z)  with  a  vertical  slope  (Figure  lb).  In 
other  words,  the  derivative  of  /  at  z  with  respect  to  arc  length  in  the  direction  given  by  6^ 
is  pure  imaginary.   Now  set 

u(,)  =  Re/(z).  (1.3) 

Then  the  derivative  of  u  in  the  direction  given  by  6^  is  zero.  Of  course  the  same  is  true  if 
/'(z)  =  0.  Thus  in  cither  case  (1.3)  satisfies  (1.1)  and  (1.2),  and  wc  have  found  a  solution 
to  Problem  O. 
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Conversely,  every  solution  to  Problem  O  is  the  real  part  of  an  analytic  function  that 
maps  0  onto  a  domain  with  straight  sides.  The  following  theorem  is  proved  in  the 
Appendix. 

THEOREM  1.  A  function  u{z)  defined  m  HU  (Ji^^i  'f  <^  solution  to  Problem  O  if  and 
only  if  u(z)  =  Ref{z),  where  f  is  an  analytic  function  in  CtU  (Ji^i  that  maps  each  interval  T^ 
onto  a  subset  of  a  straight  line  in  C  oriented  at  the  angle  9;  counterclockwise  from  the  real 
axis. 

To  illustrate  the  conformal  mapping  construction  we  will  now  look  at  an  example, 
which  appears  again  in  Section  3.  Consider  the  oblique  derivative  problem  shown  in  Fig- 
ure 2a.   The  domain  fl  is  a  rectangle  bounded  by  four  sides  Fj r4.   On  Tj  and  Tj  we 

have  Neumann  conditions,  and  on  Tj  and  V^,  oblique  derivative  conditions  at  45°  from  the 
normal  in  the  the  indicated  directions.  To  apply  the  conformal  mapping  method,  we  must 
find  a  polygonal  image  domain  /(fi)  whose  sides  are  oriented  as  suggested  in  Figure  2b. 
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Figure  2.  Problem  O  in  a  rectangle  with  82=63=0,  64= -|-,  61=- -j-. 


It  can  be  shown  that  for  this  particular  geometry,  Problem  O  has  no  bounded  non- 
constant  solutions  that  are  continuous  throughout  fl  (see  Theorem  3  below).  Instead  we 
must  look  for  solutions  that  arc  unbounded  near  one  or  more  of  Zj,  .  .  .  ,z^.  Figure  3  illus- 
trates some  of  the  possibilities. 
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Figure  3.  Various  possible  solution  domains  f(Cl)  for  the  problem  of  Fig.  2. 
In  each  case  /(fl)  lies  to  the  left  as  the  boundary  is  traversed  in  the  direction  1-2-3- 
4-1. 


First,  Figure  3a  shows  the  domain  /(O)  for  a  simple  solution  with  /(z3)  =  °o.  Here 
each  /(r^)  is  oriented  as  in  the  first  column  of  Figure  2b.  Figure  3b  shows  a  similar 
domain  with  different  orientations  of  /(rj  and  /(Fj).  Figure  3c  shov/s  a  domain  with 
/(Z2)  instead  of  /(zj)  at  00.  (Any  or  all  of  the  vertices  of  n  may  map  to  -,  but  the  C' 
hypothesis  of  Problem  O  implies  that  no  other  points  of  O  or  dCl  may  map  there.)  Figure 
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3d  shows  a  case  with  f'(z)  =  0  at  a  point  zCF^,  so  that  the  boundary  component  /(P^) 
bends  back  on  itself.  Figures  3e  and  3f  show  domains  with  two  vertices  at  <».  Figure  3g 
shows  a  domiiin  with  all  four  vertices  at  <»,  bounded  by  four  sUts,  the  cndpoint  of  each  of 
which  is  the  image  of  a  point  as  in  Figure  3d  with  f'(z)=0.  Figure  3h  shows  a  case  in 
which  /  is  not  one-to-one:  f(Cl)  is  a  tA^'o-shected  polygon. 

Obviously  the  possibilities  of  Figure  3  could  be  multiplied  ad  infinitum.  But  if  we 
exclude  functions  with  essential  singularities  at  the  vertices,  the  set  of  all  solutions  is  a 
vector  space  of  countably  infinite  dimension  that  can  be  characterized  in  an  orderly  way 
by  means  of  formulas  related  to  the  Schwarz-Christoffel  transformation.  This  is  true  not 
just  for  the  problem  of  Figure  2  but  for  the  general  instance  of  Problem  O.  We  will  now 
present  these  formulas  and  indicate  how  they  can  be  used  to  construct  solutions  numeri- 
cally. 

The  first  step,  illustrated  in  Figure  4,  is  to  transplant  the  problem  from  n  onto  the 
upper  half  f-plane  by  a  conformal  map  <}>  with  (})(Zi  )  =  «  for  some  vertex  z^  .   It  is  easy  to 

see  that  Problem  O  is  invariant  under  this  transplantation  in  the  sense  that  (j)  induces  a 
one-to-one  correspondence  between  solutions  to  Problem  O  on  ft  and  solutions  to  Problem 
O  (with  the  same  parameters  9;^)  on  the  upper  half  plane  viewed  as  a  polygon  with  degen- 
erate vertices  at  each  <t)(zi).  It  is  well  known  that  <j)  must  have  the  Schwarz-Christoffel 
form 

rKt)  =  c,  +  c,f  Uit'-Q-^^dt'  (1.4) 

for  some  constants  Ci,C2^C,  T^(.R,  and  prevertices  fi  =  4>(2i)  on  the  real  axis  [14].  Here 
P^iT  is  the  external  angle  of  8ft  at  z^,  with  -l^p^^S.  From  a  practical  point  of  view  the 
difficulty  is  that  the  values  {t^}  are  not  known  a  priori.  Determining  <}>  numerically  there- 
fore requires  an  iteration  to  find  these  values.  An  efficient  method  for  this  has  recently 
been  developed  by  the  first  author  and  implemented  in  a  Fortran  package  called  SCPACK 
[16,17].  All  of  the  computations  of  the  next  two  sections  were  performed  with  SCPACK. 
A  machine-readable  copy  of  this  package  can  be  obtained  by  writing  to  Trcfethen. 

The  second  step,  also  illustrated  in  Figure  4,  is  to  map  the  upper  half  f-plane  onto 
/(ft)  by  a  function  »j;(f).  Essentially  this  is  just  another  Schwarz-Christoffel  map  involving 
the  same  prevertices  f^  and  new  angle  parameters  3^,  but  there  are  several  new  features. 
First,  \\i  may  have  branch  points  in  the  upper  half  plane;  to  handle  these  we  must  permit 


Figure  4.   Construction  of  /  as  a  composition  of  two  Schwarz-Christoffel  maps. 

additional  factors  (t'-s^)(t'-Sj^)  in  (1.4)  [4].  Second,  additional  factors  t'-j^  with  j^  ^  R 
will  have  to  be  included  whenever  there  arc  points  along  some  F^  with  /'(z)  =  0,  as  in  Fig- 
ures 3d  and  3g.  Both  of  these  types  of  additional  factors  can  be  conveniently  absorbed 
into  a  single  polynomial  p(t')  with  real  coefficients.  Fmally,  since  f(Cl)  is  not  required  to 
be  embeddable  in  the  plane,  the  parameters  3^  are  no  longer  constrained  to  lie  in  [-1,3]. 
All  together,  the  general  map  »|;  with  algebraic  singularities  only  can  be  written  in  the 
form 


il;(f)  =  C3  +  e"'-"Jp(t')  n  it'-t,r-^''df,         Pi  =  ^(ei-ei_i), 


(1.5) 


i-1 


where  p  is  a  polynomial  with  real  coefficients,  each  m^  is  an  integer,  T^  is  any  point  in 
(ft  -1,°°),  C3  is  a  complex  constant,  and  the  branch  of  each  factor  in  the  product  is  taken 
so  that  its  value  at  T^  is  positive.  Note  that  rather  than  introduce  a  complex  constant  C4 
in  front  of  the  integral  in  analogy  with  (1.4),  we  have  prescribed  the  phase  of  this  con- 
stant explicitly  and  absorbed  its  modulus  in  p. 

The  following  theorem,  whose  proof  is  given  in  the  Appendix,  confirms  that  /  can 
always  be  constructed  from  (1.4)  and  (1.5). 
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THEOREM  2.  Let  ^  be  a  confomuil  map  (1.4)  of  D.  onto  the  upper  half  t-plane,  and 
let  (1>  be  the  function  (IJ)  defined  by  an  arbitrary  choice  ofC^.  T^,  p.  and  {mj  subject  to  the 
restrictions  listed  above.  Then  u  =  Ref  with  /=«{»o<t)  is  a  solution  to  Problem  O.  Con- 
versely, if  u  is  a  solution  to  Problem  O  that  satisfies  \u(z)\^cons,l\z-Zi^\~''^  for  some  M  in  an 
^^neighborhood  of  each  finite  vertex  z^,  and  | « (z)|:s const |z|'*^  near  any  vertex  z^  at  <»,  then 
u  =  Re/  with  f-\\io^  for  some  ij;  of  the  form  (1 .5). 

From  Theorem  2  one  can  obtain  as  a  corollary,  also  derived  in  the  Appendix,  a  com- 
plete description  of  all  bounded,  continuous  solutions  to  Problem  O: 

THEOREM  3.  Let  u  be  a  solution  to  Problem  O  that  is  continuous  on  O,  and  let  each 
6i  be  normalized  to  lie  in  (--j^'f^-  ^f  "  ^^  nonconstant,  then  it  attains  both  its  maximum 
and  its  minimum  value  at  vertices  z^  with 

— f-  <  9i  <  Oi_i  s  y,         i.e.    p4<0  (1.6) 

(and  possibly  elsewhere  also).  If  K  is  the  number  of  vertices  at  which  (1.6)  is  satisfied,  then 
the  set  of  solutions  to  Problem  O  in  C(fl)  is  the  vector  space  of  dimension  max{l,Ar}  obtained 
from  (1.5)  by  taking  m,^  =  -l  at  vertices  of  type  (1.6)  and  m,^  =  Q  at  other  vertices,  and  letting 
p  be  an  arbitrary  real  polynomial  of  degree  at  most  K—2. 

These  theorems  make  it  clear  that  the  numerical  construction  of  solutions  to  Problem 
O  should  be  manageable.  As  a  practical  observation,  note  that  once  ^  has  been  found, 
the  prcvertices  {f^}  for  ^  are  known.  Thus  in  situations  with  deg/7=0  and  m^^O  in  (1.5), 
there  is  no  paran>eter  problem  to  be  solved  for  i{>.  Our  Brownian  motion  application  in 
Section  3  is  essentially  of  this  type.  In  such  cases  the  lengths  of  the  sides  of  the  image 
domain  /(H)  cannot  be  chosen  arbitrarily,  but  arc  constrained  implicitly  by  H  and  {9^}.  In 
Figure  3  we  have  acknowledged  this  situation  by  making  each  domain  asymmetrical  even 
when  a  more  attractive  symmetrical  shape  suggests  itself  (e.g.,  let  the  two  horizontal  lines 
in  Figure  3a  lie  at  the  same  height).  Such  symmetrical  solutions  will  be  possible  in  that 
example  only  when  the  rectangle  of  Figure  2a  is  a  square. 

When  nontrivial  terms  p  and  m^  enter  into  the  map  i}>,  on  the  other  hand,  their 
values  are  normally  not  known  a  priori,  but  are  determined  by  extra  conditions  involving, 
for  example,  lengths  of  sides  of  f{il).  Our  Hall  effect  application  in  Section  2  is  of  this 
type.    In  these  cases  ^  is  no  longer  known  immediately,  but  it  can  be  determined  by 
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solving  a  generalized  Schwarz-Christoffel  parameter  problem.  Since  p  enters  (1.5)  linearly, 
the  generalized  parameter  problem  is  linear  and  easily  solvable.  The  idea  of  generalized 
parameter  problems  was  introduced  in  [18]. 

The  great  advantage  of  the  method  proposed  in  this  paper  from  a  numerical  point  of 
view,  besides  its  elegance,  is  that  it  handles  the  singularities  at  the  comers  of  fl  in  an 
essentially  exact  way.  We  expect  that  it  would  be  very  difficult  to  obtain  comparable  per- 
formance from  finite-difference  or  finite-clement  methods.  The  disadvantage  of  our 
method  is  that  it  is  very  special,  being  unable  to  treat  morc  general  domains,  boundary 
conditions,  or  differential  equations. 

Boundary-value  problems  for  elliptic  equations  in  non-smooth  domains  have  been 
studied  previously  by  a  number  of  mathematicians  from  the  viewpoint  of  the  theory  of 
partial  differential  equations.  For  surveys  with  many  references,  see  Grisvard  [5]  and 
Kondrat'ev  and  Olcinik  [11].  In  particular,  Grisvard  has  investigated  Poisson's  equation 
in  polygonal  domains  with  oblique  derivative  (and  other)  boundary  conditions.  Our  own 
work  was  motivated  originally  by  the  paper  [7]  by  Harrison,  Landau,  and  Shcpp. 
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2.  Application  1:  Hall  e£fect 

The  classical  Hall  effect  is  an  effect  that  arises  when  electric  current  flows  through  a 
conductor  or  semiconductor  in  the  presence  of  a  transverse  magnetic  field  [24].  The  mag- 
netic field  induces  a  sideways  force  on  each  moving  charged  particle,  and  this  results  in  an 
altered  current  distribution  characterized  by  a  compensating  voltage  difference  between 
one  side  of  the  conductor  and  the  other.  The  measurement  of  this  voltage  difference  is  of 
fundamental  importance  in  semiconductor  electronics.  In  particular,  its  sign  reveals 
whether  the  material  being  tested  carries  airrent  primarily  by  electrons  (type  n)  or  by 
"holes"  (type  p). 

To  be  more  precise,  let  fl  be  a  planar  polygonal  domain  as  in  Problem  O,  which  we 
will  think  of  as  having  been  cut  from  a  thin  material  of  uniform  electrical  conductivity. 
See  Figure  5.  On  each  boundary  arc  Tj  with  ;k  in  an  index  set  2C{1,  .  .  .  ,n},  a  fixed  con- 
stant voltage  u=u,^€.R  will  be  applied,  v/hile  the  remainder  of  the  boundary  is  insulated. 
We  assume  that  X  contains  no  adjacent  pairs  k,  k+1  except  with  Wi  =  «i+i.  Let  K^n/2, 
K^l  be  the  number  of  disjoint  components  of  Ui^^Fi,  that  is,  the  number  of  separate  vol- 
tages applied  on  dil. 


^ 

fai      u=u^ 

v"-^^ 

J 

!"/ 

u=u. 

i/ 

Figure  5.   The  Hall  effect  (Problem  O^)  in  a  rectangle. 
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Let  a  constant  magnetic  field  B  be  applied  perpendicular  to  O.  The  resulting  steady- 
state  voltage  m(z),  electric  field  £(z),  and  current  density  J(z)  are  related  by 

E  =  -Vu,  (2.1) 

J=-E  +  jy.S,  (2.2) 

where  we  have  set  physical  constants  to  1  for  convenience.    If  the  real  2-vectors  E  and  J 
are  interpreted  as  complex  scalars  and  B-\B\,  (2.2)  becomes 

E--{\-iB)J.  (2.3) 

Thus  E  and  J  form  a  constant  angle  a  =  tan"^B  at  every  point  z  €  fi: 

|y|/|£|  =  cosa,  argy-arg£-a,  (2.4) 

and  a  is  called  the  Hall  angle.  Without  loss  of  generality  we  may  assume  ia|<-y-. 

It  follows  from  (2.4)  that  the  problem  of  determining  u  has  the  form  of  Problem  O. 
On  each  side  with  it  €  S,  u  is  constant,  which  implies  that  it  satisfies  an  oblique  derivative 
boundary  condition  with  9^= y.  On  each  side  with  k  \^,  there  can  be  no  current  through 
the  boundary,  i.e.  /  is  tangent  to  I\.  By  (2.1),  (2.3)  and  (2.4),  Vu  therefore  makes  an 
angle  a  clockwise  from  F^  in  this  case,  which  implies  that  the  directional  derivative  of  u  is 
zero  at  an  angle  a  clockwise  from  the  normal  to  F^,  i.e.  6^=0.  Thus  an  oblique  derivative 
boundary  condition  is  indeed  satisfied  on  every  side  of  3fl.  But  in  addition  we  have  two 
further  conditions:  the  values  m^m^  on  F^  are  specified  for  fc  €  2,  and  any  physiczdly  mean- 
ingful solution  must  be  continuous.  We  combine  all  of  this  into  the  following  general  Hall 
effect  problem. 

Problem  O^.  Let  fl,  S,  u^,  and  a  be  given  as  described  above.  Find  a  solution  to 
Problem  O  subject  to  the  oblique  derivative  boundary  conditions  (1.2)  with 

64  =  f    for  it  e  2,  (2.5a) 

e^  =  a    forJtlS  (W\<f),  (2.5b) 

together  with  the  additional  conditions 

udCia),  (2.5c) 

u(z)  -  «i    for  all  z  U\,  /t  e  2.  (2.5d) 
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In  the  last  section  we  saw  that  in  general  Problem  O  has  many  solutions.  However, 
the  extra  conditions  (2.5c),  (2.5d)  of  Problem  O^  suffice  to  ensure  uniqueness.  In  the 
following  theorem  let  (J)  be  a  conformal  map  of  CI  onto  lmf>0  as  in  Section  1.  Let  {q,}, 
IsJtsAT,  be  the  subset  of  {tj  at  which  the  boundary  changes  from  type  *  ^ 2  to  type  kdl,, 
and  let  {r^},  l^k-<K,  be  the  subset  at  which  it  reverts  to  it  ^2.  In  the  product  below,  a 
factor  should  be  omitted  if  the  corresponding  point  q,^  or  r^  lies  at  <». 

THEOREM  4.  Problem  0^  has  a  unique  solution,  which  can  be  represented  in  the 
form  u  =  Re(4/o<}))  with 

m  =  C,  +  e"'-"  J  pit')Yl{f-q,y~^if-rj;)~^~^df,  (2.6) 

where  p  is  a  real  polynomial  of  degree  at  most  K—2. 

Proof.  First  consider  Problem  0^  with  condition  (2.5d)  omitted.  By  Theorem  3,  the 
solutions  are  precisely  the  functions  u  =  Re(iJ;o(}))  with  »ji  given  by  (2.6),  which  constitute  a 
vector  space  of  dimension  K.  Thus  what  we  have  to  show  is  that  the  addition  of  (2.5d) 
leads  to  existence  and  uniqueness.  Since  the  number  of  equations  in  (2.5d)  matches  the 
dimension  K  of  the  vector  space,  it  is  enough  to  show  just  one  of  these. 

We  prove  uniqueness.  If  Uj  and  «2  are  two  distinct  solutions,  then  their  difference 
«  =  «i-K2  is  harmonic  in  O,  continuous  on  fl,  zero  on  F^  for  A:  €2,  but  not  zero  every- 
where. By  Theorem  3,  the  maximum  and  minimum  of  u  on  fl  must  be  attained  at  ver- 
tices z^  satisfying  (1.6).  But  by  (2.5a-b),  the  only  such  vertices  have  Jt-1  €  2,  hence  by 
continuity  u(z J  =  0.  This  is  a  contradiction. 

An  alternative  proof  of  existence  can  be  based  on  Theorems  1.4.3.1,  1.4.4.1, 
4.4.2.1,  and  5.1.3.1  of  [5].  ^ 

The  idea  of  constructing  solutions  to  problem  O^  by  Schwarz-Christoffel  maps  goes 
back  to  an  impressive  paper  by  Wick  in  1954  [24]  and  has  been  pursued  more  recently  by 
Versncl  [21,22].  To  begin  with,  consider  the  simple  case  in  which  fi  is  a  rectangle  of 
height  1  and  length  L  as  in  Figure  5a.  Here  K=2,  so  p  reduces  to  a  constant  in  (2.6).  To 
determine  this  constant  and  Cj,  one  can  simply  set  up  a  Schwarz-Christoffel  integral  (2.6), 
then  shift  and  scale  the  resulting  parallelogram  so  that  Re/(ri)  =  Uj  and  Re/(r3)  =  «3. 
Once  this  is  done,  the  height  h  of  /(fl)  can  be  determined  by  evaluating  an  integral.  This 
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amounts  to  a  determination  of  the  total  current  /  from  Pj  to  Fj,  for  /  is  equal  to  the 
integral  of  the  magnitude  of  the  normal  component  of  /  along  Ti,  or  equivalently  along 
/(Fi),  since  current  is  invariant  under  conformal  transformation.  By  (2.4),  we  have  |£|"1 
and  |/|«cosa  along  /(Fj),  so  the  normal  component  of  /  has  constant  magnitude  cos'a. 
Therefore 


/  ~  Acos-a. 


(2.7) 


Figure  6.    Numerical  solution  of  the  Hall  effect  problem  of  Fig.  5  with  a.  =  -f- 
Streamlines  and  cquipotential  lines  are  plotted. 


Figure  6  shows  the  result  of  a  such  a  calculation  with  L  =  2,  a=^»  and  u^—u^  =  l. 
The  numerically  determined  height  and  current  are  A  =  0.615567  and  7  =  0.461675  (the 
domain  and  range  are  plotted  on  different  scales).  In  each  region  a  set  of  streamlines  and 
cquipotential  Unes  is  shown,  meeting  everywhere  at  an  angle  y-a-  The  grid  in  Figure 
6a  is  simply  the  conformal  image  in  ft  computed  by  SCPACK  of  the  grid  of  Figure  6b  in 
/(ft).   Table  1  records  computed  currents  I  for  several  other  combinations  of  parameters. 

If  ft  is  a  more  general  polygon  with  K^3  prescribed  voltages,  the  picture  changes 
qualitatively  because  now  (2.6)  contains  a  nontrivial  polynomial  p,  and  this  introduces 
sUts  and/or  branch  points  in  /(ft).  For  example,  consider  the  6-gon  shown  on  the  left  in 
Figure  7,  whose  vertices  arc  -l  +  i",  -1,  0,  -j-^j  2-4,  2  +  i.    In  this  problem  voltages 
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a 

Qo  150  30°  45^.        60^.        211  i£! 

0.5   2.000000  1.901557  1.624521  1.220042  .761746  .330865  0 

1.0   1.000000  .965926  .866025  .707107  .500000  .258819  0 

1.5    .666667  .650371  .600951  .516643  .393868  .225911  0 

2.0    .500000  .490657  .461675  .409822  .328193  .202461  0 

2.5    .400000  .393982  .375031  .340018  .281797  .183762  0 

3.0    .333333  .329141  .315803  .290595  .246974  .168283  0 

3.5    .285714  .282629  .272736  .253726  .219825  .155219  0 

4.0    .250000  .247634  .240006  .225161  .198055  .144040  0 


Table  1.    Numerically  determined  currents  I-hcosra  as  a  function  of  Hall  angle  a 
and  rectangle  length  L  for  the  problem  of  Fig.  5. 


Mi=0,  u^=Ut  and  u^=l  are  prescribed  on  sides  Fj,  F,,  and  V^,  while  Tj,  F^,  and  F^  are 
insulated.  By  Theorem  4  we  know  that  the  former  three  sides  must  be  mapped  to  vertical 
lines  and  the  latter  to  indined  ones,  and  that  the  polynomial  p  of  degree  1  in  (2.6)  may 
introduce  a  slit  in  the  boundary  of  the  sort  in  Rgures  3d  and  3g.  Let  us  write 
/>(r)  =  const (f-f.),  so  that  f.  marks  the  "extra  prevertex"  corresponding  to  the  end  of  the 
slit.  For  arbitrary  U,  t.  might  lie  anywhere  on  R,  but  let  us  assume  0<C/<1.  Then  there 
are  exactly  three  possible  topologies  for  /(H),  and  these  are  shown  in  Figure  7.  In  each 
case  one  of  the  sides  of  /(O)  protrudes  beyond  the  vertex  it  might  be  expected  to  stop  at, 

then  doubles  back,  so  that  /(fl)  is  not  a  6-gon  but  a  7-gon,  with  vertices  w^ Wf,  and 

w..    The  three  cases  differ  in  the  order  in  which  the  vertices  and  prevertioes  appear 

around  the  boundary: 

\ 
(a):     H/j,  W2,  w.,  w^,  W4,  w^,  Wf, 

(b):     wy,  wi,  W3,  »v-,  w^,  w^,  Wf, 

(c):       Wi,  Wn,   VVj,   W^,   W.,  W5,  Wf, 

They  also  differ  physically: 
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1 
2 


u=0 


u=l 


(a) 


(b) 


Figure  7.   Hall  effect  on  a  6-gon  -  introduction  of  a  slit  of  unknown  length  in  /(H). 


(a):    current  flows  into  r3  only, 

(b):    current  flows  both  into  and  out  of  Tj, 

(c):    current  flows  out  of  r3  only. 

Of  course  it  may  happen  that  /(H)  has  no  slit  in  it,  but  only  when  the  problem  by  chance 
falls  on  the  borderline  between  (a)  and  {b)  or  (b)  and  (c). 

If  (/  is  specified  in  Figure  7,  we  do  not  know  a  priori  which  of  the  cases  (a)-(c)  will 
occur,  nor  do  we  know  the  correct  slit  length  or  any  of  the  other  vertical  dimensions  of 
/(H).  On  the  other  hand  we  do  know  its  horizontal  (i.e.  oblique)  dimensions.  This  is 
why  the  determination  of  i[i  is  a  generalized  parameter  problem,  in  contrast  to  the  usual 
Schwarz-Christoffel  situation  in  which  aU  dimensions  of  the  image  polygon  arc  known  at 
the  start.    But  since  the  generalized  parameter  problem  is  linear,  it  is  easy  to  solve  despite 
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the  complexity  of  its  geometric  interpretation.  A  related  but  nonlinear  generalized  param- 
eter problem  is  described  in  the  paper  in  this  volume  by  Herat  and  Trefethcn  on  free- 
streamline  flows  [3]. 


1    ;    I    i    1    /    1    , 

'     /    /   J^-l—4—^~~r't- 

^^^ 

"  /-'v — 1    1    /    -^    /    /  — ' 

-^  L^—r — 1    /    /     /     /     /     / 

'J -'T — Z__  /      /      /—/      '      I       1      /— ~y — t     /-W'    ^  1 

A'A'^i^  /  '   "^   /  ! 

^m 

^^f-^-f^^^y^^^ 

X^ 

Figure  8.  Numerical  solution  of  the  Hall  effect  problem  of  Fig.  7,  a=-^. 


For  a  practical  illustration  with  K>2,  we  computed  the  numerical  solution  shown  in 
Figure  8  with  t/=-,  a=T--  ^^^  ^*^*  ^^  "^  '^  '^^'^  ^  solution  has  turned  out  to  be  of 
type  {h).  The  current  between  each  adjacent  pair  of  streamlines  is  0.05.  The  total  currents 
arc  .449016  out  of  Tj,  .166622  into  Tj,  .045959  out  of  Tj,  and  .328354  into  r^. 

The  numerical  method  we  have  described  in  this  section  is  new  even  in  the  case 
5=a=0  in  which  the  magnetic  field  is  not  present,  so  all  that  remains  is  a  mixed 
Dirichlet-Neumann  problem.  Now  m  is  the  real  part  of  a  conformal  map  /  onto  a  domain 
with  horizontal  and  vertical  sides  and  with  sides  and  slits  of  unknown  length.  To  illus- 
trate, Figure  9  shows  a  computation  for  the  same  problem  as  in  Figure  8  but  with  a=0. 
Again  the  solution  turns  out  to  be  of  type  (b).  Now  the  total  currents  arc  .452474  out  of 
Tj,  .D6224  into  Tj,  .079972  out  of  Tj,  and  .396222  into  Fj. 

To  show  that  cases  {a)  and  {c)  can  indeed  occur,  Figure  10  presents  a  3x3  array  of 
numerically  determined  boundary  polygons  /(afi)  as  in  Figures  8  and  9  for  the  problem  of 
Figure  7  with  a  =  --^,0,  y  and  t/=.l,.5,.9.  In  this  experiment  wc  have  evidently 
obtained  case  {a)  with  t/  =  .l,  (h)  with  t/=.5,  and  {c)  with  t/  =  .9. 
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1 — 1- 

Figure  9.  Numerical  solution  of  the  Hall  effect  problem  of  Fig.  7,  a  =  0. 


Most  of  the  computer  time  in  any  SCPACK  computation  is  ultimately  spent  in  the 
calculation  of  complex  logarithms  needed  to  evaluate  the  powers  of  (f'-f^)  in  (1.4)  and 
(1.5).  As  a  result,  one  can  conveniently  estimate  the  work  done  by  SCPACK  in  a 
machine-independent  way  by  counting  logarithms.  In  practice  the  most  time-consuming 
parts  of  the  computation  arc  the  solution  of  nonlinear  parameter  problems,  where  the  cost 
increases  roughly  in  proportion  to  n^,  and  the  construction  of  plots,  where  mapping  of 
many  points  is  required,  y^jproximatc  logarithm  coimts  for  the  computations  of  this  sec- 
tion are  90,000  for  Figure  6,  80,000  for  Table  1,  180,000  each  for  Figures  8  and  9,  and 
20,000  for  Figure  10. 

In  their  papers  on  the  Hall  effect,  Wick  and  Versnel  consider  the  situation  in  which 
two  edges  are  held  at  fixed  voltages  while  several  others  are  passive  conductors  whose 
voltage  is  free  to  vary,  but  across  which  no  net  current  may  flow  [21,22,24].  This  leads 
to  a  somewhat  different  linear  generalized  parameter  problem  in  which  some  of  the  equa- 
tions enforce  the  condition  that  the  two  sides  of  a  slit  must  have  equal  (but  unknown) 
length.  Wick  and  Versnel  confine  themselves  to  symmetrical  domains  for  which  the  map 
<})  can  be  obtained  more  or  less  analytically. 
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Figure  10.   Image  polygons  /(afi)  for  the  Hall  effect  problem  of  Fig.  7  for  various  a 
and  U. 
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3.  Application  2:  reflected  Brownian  motion 

The  mathematical  concept  of  two-dimensional  Brownian  motion  is  a  foimalization  of 
the  idea  of  random  motion  of  a  particle  in  the  plane  [10].  Under  suitable  restrictions, 
Brownian  motion  confined  to  a  domain  can  be  described  by  a  probability  density  p(z,t) 
that  diffuses  with  t  according  to  the  heat  equation,  and  in  the  limit  f-oo  converges  to  a 
steady-state  probability  density  u{z)  that  satisfies  Laplace's  equation.  Of  course  the  nature 
of  u  will  depend  on  boundary  conditions,  which  describe  the  behavior  of  the  particle  when 
it  hits  the  boundary.  The  connection  of  this  topic  with  the  present  paper  is  that  in  several 
problems  of  practical  interest,  the  appropriate  boundary  conditions  for  the  steady-state 
density  involve  oblique  reflection,  which  leads  to  instances  of  Problem  O.  In  particular,  in 
queueing  theory,  Brownian  motion  with  oblique  reflection  at  the  boundary  of  a  polygonal 
domain  arises  as  a  diffusion  approximation  to  the  queue-length  process  of  a  tandem  queue 
under  conditions  of  heavy  traffic  (approximate  equality  of  average  arrival  and  service 
rates)  [6].  The  same  problem  arises  under  similar  conditions  in  storage  theory  in  modeling 
the  contents  process  of  a  two-buffer  production  or  storage  network  [2,23]. 

To  be  more  precise,  let  fl  and  Bj,  .  .  .  ,9„  be  defined  as  in  Problem  O,  with  |6il<Y 
for  all  it.  The  goal  is  to  construct  the  steady-state  density  u  of  a  strong  Markov  process  Z 
with  continuous  sample  paths  that  (a)  behaves  like  Brownian  motion  in  the  interior  of  ft, 
and  (b)  is  reflected  instantaneously  on  each  side  T^  at  the  angle  -9^  from  the  normal.  It 
is  known  that  such  a  process  can  always  be  constructed  with  absorption  at  the  vertices  of 
ft.  However,  if  ft  is  a  bounded  domain  with  external  angle  parameters  p^  in  the  range 
(-1,1),  and  if  the  reflection  does  not  provide  too  great  a  "push"  towards  any  finite  ver- 
tex, then  it  has  also  been  shown  that  there  is  also  a  unique  such  process  that  spends  zero 
time  (in  the  sense  of  Lebesgue  measure)  at  the  vertices.  A  necessary  and  sufficient  condi- 
tion for  this  is 

Pi  <  1  -    ^~  ^"^      for  all  k.  (3.1) 

The  formal  definition  of  the  reflected  Brownian  motion  process  and  the  derivation  of 
these  results  follow  from  the  work  of  Varadhan  and  Williams  [19]  by  localization,  as 
explained  by  Harrison,  Landau,  and  Shepp  in  [7]. 

For  example,  consider  again  Figure  2a.  Here  62=63=0,  Qi=-6i==-^,  and  Pi  =  Y  for 
k-1,2,3,4.  The  reflected  Brownian  motion  for  this  geometry  can  be  interpreted  as  the  dif- 
fusion limit  of  a  system  of  two  finite-capacity  queues  Qi  and  Q2  in  parallel,  where  the  x 
and  y  positions  in  the  rectangle  represent  the  normalized  numbers  of  customers  waiting  in 
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Cj  and  Q2,  respectively,  and  the  queues  are  coupled  as  follows:  if  Q^  becomes  empty,  it 
takes  customers  from  Q2,  while  if  Q2  becomes  full,  it  transfers  customers  to  Ci-  Since 
(3.1)  is  satisfied,  there  is  a  unique  process  Z  that  spends  zero  time  at  the  vertices,  behaves 
like  Brownian  motion  inside  the  rectangle,  and  is  reflected  instantaneously  in  the  normal 
direction  on  Tj  and  Fj  and  at  angles  -^  to  the  right  on  r4  and  —■  downwards  on  V^.  Note 
that  in  keeping  with  the  minus  sign  in  condition  (b)  above,  the  directions  of  reflection  on 
each  side  are  not  the  same  as  the  directions  in  which  the  derivative  is  zero,  but  are  equal 
to  these  directions  flipped  about  the  respective  normals. 

Our  purpose  is  to  compute  the  steady-state  density  u  for  reflected  Brownian  motion 
on  a  bounded  polygon  Ci  by  taking  the  real  part  of  a  conformal  map  /  of  the  kind 
described  in  Theorems  1  and  2.  As  usual,  the  statement  of  Problem  O  does  not  by  itself 
include  enough  conditions  to  characterize  u.  The  additional  conditions  needed  turn  out  to 
be  quite  different  from  those  that  made  sense  for  the  Hall  effect  problem.  First,  since  u  is 
a  probability  density  it  must  be  nonncgative  -  f(il)  must  lie  in  the  right  half  plane. 
Second,  the  total  probabihty  J  u  must  be  finite  and  equal  to  1  --  this  determines  the  scale 
for/(n).  Finally,  there  is  a  third  restriction  that  is  less  obvious.  This  is  a  global  "comer 
condition"  that  amounts  to  a  statement  that  particles  may  not  disappear  at  one  vertex  and 
reappear  at  another: 

Comer  condition,  ffa  and  b  are  two  points  on  any  two  sides  Vj  and  Fj,  respectively, 
with  Q{a)  =  Qj  and  Q(b)  =  Qi^,  and  A  is  a  differentiable  curve  in  CL  that  joins  a  to  b,  then 

-f—(.z)da=  «(i)tane(fe)-u(a)tan9(a),  (3.2) 

A  ^^ 

where  -^  denotes  the  derivative  in  the  normal  direction  to  A  that  lies  to  the  left  as  A  is 

an 

traversed  fi^om  a  to  b. 

Expressed  in  this  form  the  comer  condition  may  seem  obscure,  but  there  is  a  simple 
equivalent  statement  in  terms  of  the  analytic  function  /:  it  must  be  possible  to  choose  /  so 
that  /(n)  is  a  radial  domain.  By  this  we  mean  that  all  of  the  boundary  arcs  /(F^)  are  sub- 
sets of  straight  lines  passing  through  the  origin,  hence  in  this  particular  application, 

arg/(z)  ^  6i  (mod  it)       for  aU  z  €  F^.  (3.3) 

This  observation  was  first  made  by  Newell   [15].     To  sec  that   (3.2)  and   (3.3)  arc 
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cqnivalcnt,  note  that  if  /=u  +  jv,  then  by  the  Cauchy-Riemann  conditions,  -j^  =  -j^- 

Therefore  the  left  side  of  (3.2)  is  equal  to  v(fe)-v(a),  so  (3.2)  anoounts  to  the  statement 
that  v(z)  =  M(z)tane(z)  +  C  on  each  P^,  for  sonie  fixed  constant  C.  By  choosing  v  so  that 
this  constant  is  zero,  (3.2)  is  reduced  to  (3.3). 

The  determination  of  conformal  maps  onto  radial  domains  is  easier  than  the  general 
Schwarz-Christoffel  problem,  since  arg/  rather  than  just  arg/'  is  known  on  each  side  T^. 
By  repeating  the  proof  of  Theorem  2  with  appropriate  changes,  one  can  show  that  assum- 
ing u  grows  at  most  algebraically  at  each  vertex,  /  can  be  written  in  the  form  /=i{»o<j) 
with 

m  =  A"p(o  n  it-t,r~^'  (3.4) 

i-l 

in  analogy  to  (1.5).  Thus  no  integration  is  required  for  ij;,  which  is  of  course  a  consider- 
able help  in  numerical  computations.  We  will  see  further  in  Theorem  5  below  that  if  «  is 
nonnegative,  then  p  and  {m^}  in  (3.4)  can  be  omitted. 

Returning  to  the  example  of  Figure  2,  recall  that  many  solutions  to  Problem  O  were 
depicted  in  Figure  3.  In  the  present  situation  all  of  these  except  Figure  3f  are  excluded  by 
the  nonncgativity  condition  on  m.  Figure  3f,  in  turn,  is  excluded  by  the  condition  that 
f(Cl)  must  be  radial.  To  get  a  radial  domain,  the  slit  ^2-^3-^4  in  Figure  3f  must  be 
lowered  to  the  real  axis.  Unless  n  is  a  square,  w^  will  have  to  move  away  from  the  end 
of  the  slit  to  make  this  possible.  Thus  the  proper  solution  to  the  reflected  Brownian 
motion  problem  of  Figure  2  is  given  by  the  conformal  map  indicated  in  Figure  11.  The 
new  symbols  z.  and  w.  denote  the  end  of  the  sMt  and  its  preimage  on  dCt. 

We  are  now  ready  to  formulate  the  variant  of  Problem  O  motivated  by  the  reflected 
Brownian  motion  problem.  Note  that  in  contrast  to  Problem  0^.  here  we  do  not  require 
u  e  C(ft),  and  indeed  often  u  will  be  mfinite  at  a  vertex. 

Problem  Og.  Let  il  and  6^  be  given  as  usual  with  |6i|<Y /or  each  k.  Find  a  solution 
to  Problem  O  subject  to  the  additional  conditions 

(0      u{z)^Q     for  all  z  i  ^V^  \J  j: i,; 

(ii)     u^L\n)     with     /^"  =  i; 
(Hi)      the  comer  condition  (3.2)  is  satisfied. 
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Figure  11.   Reflected  Brownian  motion  (Problem  Ob)  in  a  rectangle. 

Like  Problem  O^ ,  Problem  Og  is  well  posed.   As  usual,  in  the  following  theorem  cj) 
is  a  fixed  conformal  map  (1.4)  of  £t  onto  lmr>0. 

THEOREM  5.    Problem  O^  has  a  solution  if  and  only  if  (3.1)  holds  for  each  k.    The 
solution  is  unique,  and  can  be  represented  in  the  form  K=Re(»l;oij))  with 


HO  =  ue  *-"  n  0-O"^       Pi  =  ^(Q.-h-i),  V  e  (0,00). 


(3.5) 


Proof  Suppose  u  is  a  solution  to  Problem  Og.  By  Theorem  1,  u  can  be  written  as 
tt  =  Re/=Re(iJ;o(}))  for  some  function  i|;  analytic  in  Imr&O  except  possibly  at  the  points  {rj. 
By  conditions  (i)  and  (Hi)   (=(3.3)),  ^  must  map  an  interval  [fi-£,/i)  into  the  ray 

[0,oo)e'*"S  an  interval  ('tJi+e]  into  the  ray  [0,'^)e  *,  and  the  half-plane  lmf>0  into  the 
half -plane  [argwl^  y.   Therefore  the  function 

g,{t)  =  Ht)(t-t,f'' 

maps  both  of  these  intervals  into  [O.oo)^  '.  If  ^k-\-^k.f  ^^^^  values  assumed  by  g^  in  Im/>0 
satisfy  -Ysarg>v< Y  +  ei-e^-i,  while  if  Qi.^^^%;,,  they  satisfy  -^  +  ei-0i_i<argH'<-5^. 
Since  |9;,_,|<-|^  by  assumption,  i.e.,  |8i_i|==^~S  with  5>0,  wc  have  in  cither  case 


|arg^i(r)-(ei±Ti)|s5>0 
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for  lraf>0.  In  other  words,  g^  takes  no  values  in  the  sector  of  half-angle  8>0  about  the 
ray  -[0,t»)e'  ». 

Now  gi  can  be  extended  to  lraf<0  by  reflection  in  either  [ti,-e,t^)  or  (ft.ri+e].   Since 

both  of  these  intervals  map  into  the  same  ray  [0,«>)e  ',  the  extension  is  independent  of 
which  interval  is  chosen.  Thus  g^  is  a  single-valued  analytic  function  in  a  punctured 
neighborhood  of  f^,  and  must  have  an  essential  singularity,  a  pole,  or  a  removable  singu- 
larity at  f^.  But  by  the  geometry  of  the  reflection  process  it  is  clear  that  the  extended 
function  g^  also  takes  no  values  in  the  sector  described  above.  Therefore  f^  can  only  be  a 
removable  singularity. 

These  arguments  imply  that 

m  n  ('-'/' 

i-i 

can  be  extended  by  reflection  in  R  to  a  function  that  is  analytic  for  all  r  €  C.  Moreover, 
by  reasoning  as  above  it  follows  that  it  is  also  analytic  at  ?  =  =».   Therefore  it  is  a  constant, 

and  since  its  argument  is  e   '  ,  this  establishes  the  representation  (3.5). 

At  this  point  we  have  shown  that  any  solution  u  to  Problem  Oq  has  the  form 
claimed.  Conversely,  given  an  instance  of  Problem  Og,  let  «  be  constructed  as  Re((j;o<}>) 
from  (3.5).  It  is  easy  to  sec  that  conditions  (i)  and  (Hi)  arc  satisfied.  Turning  to  (u),  it  is 
clear  that  u  will  be  integrable  on  ft  if  and  only  if  it  is  integrablc  near  each  2^.  If  z^  is  fin- 
ite, /  has  a  singularity  of  type 


\z-z,\    -'^-^'^ 

there,  and  u  is  integrable  if  and  only  if  the  exponent  is  greater  than  -2,  i.e.,  if  and  only  if 
(3.1)  holds.  A  similar  consideration  gives  (3.1)  also  for  Zi  =  <».  This  establishes  the 
existence  claim  of  the  theorem.  Since  the  parameter  U  must  be  chosen  to  ensure  J"  «  =  1, 
the  solution  is  unique. 

By  a  slight  modification  of  the  results  of  [7]  and  [25],  it  follows  that  the  solution  u  of 
Problem  Og  given  by  Theorem  5  is  the  unique  stationary  density  for  the  Markov  process 
Z,  provided  that  ft  is  bounded  with  (3^  €  (-1,1)  for  each  k.  As  conjectured  in  [7]  but  nnt 
yet  proved,  the  same  result  probably  holds  for  unbounded  polygons  with  Pie(-l,v)- 
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With  ?i^\  there  is  no  chance  of  satisfying  (3.1). 

Wc  will  now  illustrate  the  Schwarz-Christoffel  solution  of  Problem  Og  by  giving 
results  computed  with  SCPACK  for  the  problem  of  Figure  11.  (Actually  this  geometry  is 
simple  enough  that  it  could  have  been  treated  semi-analytically  via  elliptic  functions.)  Let 
fi  have  height  1  and  length  L.  In  /(H),  any  vertical  line  Re/=  const  is  a  line  of  constant 
probability  density.  Figure  12  shows  /(H)  together  with  a  system  of  such  lines 
M  =  .1,  .2,  •  ••  ,  2  for  the  case  L  =  2.  By  mapping  these  lines  back  to  O,  we  obtain  curves 
of  constant  probability  density  in  the  problem  domain  which  reveal  quite  a  bit  about  the 
steady-state  behavior  of  the  reflected  Brownian  motion.  Figure  13  shows  results  of  this 
kind  for  1  =  1,2,3.  Beginning  at  the  upper-left  comer,  we  have  m=0  at  Zj,  then  the  curves 
«  =  .l,  .2,  .  .  .  ,  2  in  succession.  Note  that  u  has  poles  at  both  z,  and  z^,  but  that  as  L 
increases,  the  former  rapidly  weakens  in  strength.  The  explanation  for  this  is  that  the 
rightward  motion  imparted  to  the  process  by  the  upper  side  makes  the  probability  density 
decay  exponentially  with  distance  from  the  right  side.  At  L=3  the  right-wall  behavior  is 
already  very  close  to  what  it  would  be  in  the  limit  L=oo,  a  case  investigated  earlier  by 
Harrison  and  Shepp  in  [9].  In  this  limit  w.  will  approach  0,  and  the  lower  half  of  /(ft) 
will  correspond  to  a  small  portion  of  ft  near  Z2. 

Similar  plots  showing  the  behavior  of  u  near  the  comers  have  been  constructed  previ- 
ously by  Newell  by  hand.  See  Figures  IV- 13  to  IV- 15  of  [15]. 

Table  2  lists  computed  values  of 


/,  =  fu,      Jt=l,2,3,4  (3.6) 

and 

<j:>  =  fxu(z),  <y>  =  Jyu^z)  (3.7) 

n  n 

for  various  lengths  L,  together  with  corresponding  quantities  w,  and  ^3.  Hfcre  the 
integrals  along  F^  are  taken  with  respect  to  arc  length  and  those  over  ft  arc  taken  with 
respect  to  Lebcsgue  measure  in  R^.  The  number  /^  can  be  interpreted  as  a  measure  of  the 
rate  at  which  the  Markov  process  Z  hits  side  F^  in  the  steady  state,  or  alternatively,  as  a 
measure  of  the  average  amount  of  "pushing"  that  must  be  exerted  on  T^  to  keep  the 
Brownian  motion  in  ft.  <j:>  and  <y>  represent  the  long-run  average  x  and  y  positions. 
In  practice  these  quantities  are  most  conveniently  computed  as  onc-dimcnsional  integrals 
by  means  of  the  identities 
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Figure  12.  Image  domain  /(O)  for  the  reflected  Brownian  motion  problem  of  Fig. 
11  with  Z,  =  2.  The  vertical  lines  represent  probabihty  densities  u  =  .1,  .2,  •■  •  ,  2,  and 
the  arrow  indicates  w,. 
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Figure  13.  Nunvirical  solution  of  the  reflected  Brownian  motion  problem  of  Fig.  11 
with  L  =  1,2,3.  The  curves  are  lines  of  constant  probabihty  density  u  =  .1,  .2,  •  •  •  ,  2 
obtained  by  applying  /"'  to  Fig.  12  or  its  analogs  for  Z,  =  1,3. 
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<x>  =  lL'^I,-Lfux\  <y>  =  i/4+iJ«r-  (3.8) 

These  relations  were  obtained  by  substituting  g=x^  and  ^=>'  in  the  following  equation 
derived  by  applying  Green's  theorem  to  smooth  domains  approximating  CI  and  taking  lim- 
its, using  property  (3.2)  of  u  and  the  behavior  of  u  near  comers: 

I«  A^  +  i  /   uD,8  =  0       for  all  g  e  C^U),  (3.9) 

where 

Di,g  =  g„-tanQi^g„,     fc=l 4. 

Similarly,  the  magnitude  of  the  constant  U  in  (3.5)  needed  to  accomplish  the  normaliza- 
tion (ji)  can  be  determined  by  computing  f  u  as  the  one-dimensional  integral 

Su=h  +  fyu  (3.10) 

derived  by  substituting  g=y^  in  (3.9).  For  numerical  evaluation  of  (3.8)  and  (3.10)  a 
natural  method  is  Gauss-Jacobi  quadrature,  since  the  integrands  are  analytic  except  for 
algebraic  singularities  of  known  type  at  the  endpoints.  This  approach  conforms  nicely  to 
the  spirit  of  SCPACK,  where  all  Schwarz-Christoffel  integrals  are  already  computed  by 
compound  Gauss-Jacobi  quadrature  using  a  well-known  program  by  Golub  and  Welsch 
referenced  in  [16]. 

By  substituting  g=x  and  g=y  in  (3.9),  one  obtains  the  interesting  balance  relations 

I^  =  I,  +  I,  =  I^  (3.11) 

(compare  [15],  eq.  in.2.3f).  The  left  and  right  equalities  here  mirror  the  fact  that  in  the 
steady  state,  there  can  be  no  net  force  of  the  sides  on  the  Brownian  motion  in  the  vertical 
and  horizontal  directions,  respectively.  These  identities  are  verified  numerically  by  the 
data  in  Table  2.  The  final  line  in  the  table  lists  the  limiting  value  of  each  column  as  L-oo, 
where  H  becomes  the  half-strip  considered  in  [9].  The  limit  for  ^3  has  been  calculated 
previously  in  eq.  (1.2)  of  [9], 

lira  W3  =  :n_2TL  ~  599070. 

r(i)2 


The  limit 


lira  L-<x>  ~  1.242454 


26- 


W3        I^      12=13        I4 


0.5  1.693474  2.912616  1.871418  2.412890  0.541472  0.258585  0.380018 

1.0  1.138859  1.138859  0.805295  1.610589  0.805295  0.551506  0.448494 

1.5  0.686627  0.831998  0.446669  1.340225  0.893557  0.879534  0.471624 

2.0  0.423368  0.728154  0.270736  1.206445  0.935709  1.239964  0.482830 

2.5  0.268692  0.678353  0.171214  1.130587  0.959373  1.628342  0.489146 

3.0  0.174154  0.649951  0.110891  1.084582  0.973691  2.040075  0.492970 

3.5  0.114465  0.632412  0.072873  1.055585  0.982712  2.471022  0.495381 

4.0  0.075922  0.621165  0.048334  1.036868  0.988534  2.917572  0.496936 


0.599070       Oil     L-1. 242454 


Table  2.   Numerically  determined  quantities  w,,  w^,  4,  <j:>,  <y>  as  a  function  of 
rectangle  length  L. 

was  obtained  numerically  with  SCPACX  by  treating  the  oblique  derivative  problem  for  a 
semi-infmite  strip  directly.  Each  column  in  Table  2  apparentiy  approaches  its  limit 
geometrically,  with  an  error  asymptotically  proportional  to  exp(-T:L/4). 

As  in  the  last  section,  we  can  measure  the  bulk  of  the  effort  required  for  these  com- 
putations in  a  machine-independent  way  by  counting  complex  logarithms.  Roughly,  Fig- 
ures 12  and  13  required  70,000  logs  and  Table  2  required  180,000. 


1 
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Appendix:  basic  theorems  on  the  oblique  derivative  problem 

In  this  appendix  wc  give  three  proofs  that  were  omitted  from  Section  1  in  the  interest 
of  readability.  The  first  argument  confirms  that  all  solutions  to  Problem  O  can  be 
obtained  by  conformal  mapping. 

Proof  of  Theorem  1.  The  proof  of  the  "if'  part  is  straightforward,  and  was  sketched 
in  the  discussion  in  Section  1.  To  prove  the  "only  if'  part,  let  m  be  a  solution  to  Problem 
O.   Since  m  e  C-(n)nc'(nu  (Ji'^i).  the  function 

g{z)  =  uXz)-iu,{z)  (A.1) 

belongs  to  C''(n)nC^'(n U  (J^ri),  and  since  u  is  harmonic  in  H,  g  satisfies  the  Cauchy- 
Riemann  equations  and  is  therefore  analytic  in  Cl.  Suppose  F^  is  inclined  at  an  angle  Tq^ 
counterclockwise  from  the  real  axis.  By  an  elementary  calculation  based  on  (1.2)  and 
(A.1),  one  can  verify  that 

arg^C^)  -  ^^-Tii  (mod-n)      for  all  zd  F^.  (A.2) 

Therefore  g  maps  each  straight  boundary  segment  F^  into  a  subset  of  the  straight  line  in 
the  complex  plane  passing  at  angle  e^-Ti^  through  the  origin.    By  the  reflection  principle, 
it  follows  that  g  can  be  extended  analytically  across  F^.   Thus  g  is  analytic  in  HU  U/t^^i- 
Let  z-  G  n  be  arbitrary,  and  consider  the  new  function 

r 

f(z)  =  u(z,)  +  fg(z')dz',  (A.3) 

I. 

where  the  path  of  integration  is  any  rcctifiable  arc  in  flu  Ui^i  ^^^  pi"s  z-  to  z.  Since  g 
is  analytic  and  fl  is  simply  connected,  /(z)  is  single-valued  and  analytic  throughout 
fiU  U^F^.    Moreover  for  z  6  HU  Ui^i'  ^^  ^'^^ 

u(z)  -  u{z.)  +  fu^dx  +  u^.dy 
I. 

r 

=  u(z.)  +  jRe{{u,-iuXdx+idy)}  =  Re/(z). 
Now  for  any  two  points  z,,z2  d  l\,  (A.3)  implies 

■2 
/(Z;)-/(Z.)     -    Jgiz)dZ. 
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Sincc  dz  has  constant  argunncnt  Ti^^CinodTr)  along  \\,  it  follows  from  (A. 2)  that /(zjj-/(z,) 

has  argument  ^X^oAti).    Therefore  /  maps  l'^  int(;  a  straight  line  oriented  at  angle  0^ 

counterclockwise  from  the  real  axis,  as  claimed.  „ 

D 

The  next  argument  shows  that  /  can  tx;  written  in  Schwarz-Christciffel  form. 

Proof  of  Theorem  2.  Suppose  first  that  ij;  is  a  map  of  the  form  05).  liy  construc- 
tion, the  integrand  of  (1.5)  (including  the  factor  in  front)  has  argument  6^     ;  (modTr)  at 

T^  and  therefore  on  4)(rj  _.).  Following  the  changes  in  argument  of  each  factor  O'-fj  * 
as  t'  moves  leftward  along  the  real  axis,  we  sec  that  the  argument  is  0^,  (modTr)  on  each 
interval  i^{\\).  Therefore  i>  is  an  analytic  function  in  lrar>0  that  maps  each  W\)  into  a 
subset  of  a  straight  line  oriented  at  the  angle  9^  from  the  real  axis.  By  Theorem  1 ,  it  fol- 
lows that  u  =  Refi^oi}))    is  a  solution  to  Problem  O. 

Conversely,  let  u  be  a  solution  to  Problem  O  satisfying  the  given  bounds  at  each  ver- 
tex. A.ssume  without  loss  of  generality  that  we  have  already  transplanted  the  problem  v) 
that  fi  is  the  upper  half  /-plane;  since  the  Schwarz-Christoffcl  map  <J)  has  only  algebraic 
singularities  at  the  comers  of  O,  this  will  not  interfere  with  the  gnrwth  condition.  For 
simplicity  assume  further  (a  minor  modification  of  (1.4))  that  ti^-^{z^)  is  finite  for  each  k. 
Let  ^  be  the  analytic  function  with  w^ReiJ/  provided  by  Theorem  1.  CXir  goal  is  to  show 
that  4(  can  be  written  in  the  form  (1.5). 

The  first  task  is  to  show  that  the  algebraic  growth  bounds  on  u  imply  that  »J/'  satisfies 
similar  bounds.  To  this  end,  assunx;  /;  =  0  for  convenience  and  suppose  |«(f)|5:const|r|~''', 
A/>0.    For  t=r+is  near  0  with  5>0,  we  can  compute  i|)'  as  the  integral 

obtained  by  differentiating  the  Schwarz  formula  [1,  p.  167].   Therefore  |^^'|  is  bounded  by 

|ii/'(Oi  ^const|r|-'^5-^.  (A.4) 

In  any  sector  bounded  away  from  the  real  axis,  say  for  jay|r|,  s~'  is  equivalent  to  \l\~\ 
and  (A.4)  implies 

|.|/'(/)|r;  constlrl--'^-'.  (.\..5; 

On  the  other  hand  consider  /  v/ith  *<-7|r|.    Integrating  i|i'  up  to  t  from  sonrx:  fixed  point, 


90 


wc  obtain  from  (A.  4)  the  estimate 


log|j|-log|f| 


(A.6) 


|»i;(0|  s  const|r|~'*' 

for  |r|  sufficiently  small.  Note  that|log|j|-log|r||>  log2  for  |j|<y|r|.  (The  term  -log|r| 
comes  from  the  assumption  A/>0,  which  ensures  that  the  integral  of  iV'  up  to  the  point 
r  +  i|/|  has  order  {tl'"^,  not  |r|"'*^log|/|.)  Now  by  Theorem  1,  »Ji  maps  a  portion  of  the  real 
line  near  t  into  a  straight  line.  In  particular,  i|»  can  be  extended  by  symmetry  to  a  func- 
tion satisfying  (A.6)  on  and  inside  the  circle  {t:  |T-f|  =  y|r|},  and  i|;'  can  then  be  obtained 
by  a  Cauchy  integral  over  this  circle.  If  wc  estimate  this  integral  from  (A.6),  the  log  fac- 
tors drop  out  in  the  integration,  and  we  arc  left  with  (A.5)  for  all  r,  as  desired. 

The  remainder  of  the  proof  is  a  modification  of  the  usual  Schwarz-Christoffel  argu- 
ment based  on  the  symmetry  principle  [4,14].  Let  r^  be  a  prcvertex.  By  Theorem  1,  i|»' 
mapjs  an  interval  [fi-t,ri)  into  the  line  through  the  origin  oriented  at  the  angle  6i,_,  from 
the  real  axis,  and  an  interval  ('i,fi+t]  into  the  line  through  the  origin  at  angle  0^,.  There- 
fore for  any  integer  m^,  the  function  >|;'(f)(r-fj^*  "'  maps  [ft-e.r^)  and  (f^.f^+t]  into  the 
single  line  e  'R.  Since  i(»'  has  at  most  algebraic  growth  at  f^,  we  can  choose  m^  so  that  this 
function  is  bounded  and  continuous  near  r^.  By  the  reflection  principle,  it  follows  that  it 
has  an  analytic  extension  to  a  disk  around  z^. 

These  observations  and  choices  of  m^  imply  that  the  function 


pit)  =  e~  '-"'^'(0  n(r-r>"'"' 


is  analytic  throughout  lraf>0  and  maps  R  into  itself.  Moreover  at  r  =  oo,  ij;  is  analytic  by 
Theorem  1  and  our  choice  of  4),  so  p  has  at  most  a  pole  (since  2(3^  =  0).  By  reflection 
across  R,  p  can  be  extended  to  an  entire  function.  Thus  p  must  be  a  polynomial  with  real 
coefficients.    Eq.  (1.5)  now  follows  by  integration. 

The  final  argument  shows  that  the  bounded,  continuous  solutions  to  Problem  O 
make  up  a  vector  space  of  finite  dimension. 

Proof  of  Theorem  3.    By  Theorem  2,  any  continuous  (hence  bounded)  solution  to 
Prnbirm  O  can  be  written  w  =  Refi];o(}))  with  ij;  given  by  fl.5).    At  any  finite  prcvertex  r, 

the  integrand  of  (1.5)  has  a  singularity  of  the  form  (r-r)^'  with  71^/7)^-3;  +  /;,,  where 
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n^^O  is  the  number  of  zeros  of  p  at  r^,  if  any,  and  u  will  be  continuous  at  r^  if  and  only  if 
"yi>-l.  Therefore  we  can  obtain  all  continuous  solutions  by  picking  m^  so  that 
'"i~Pi^(~1.0]>  ^^^  is,  by  taking  mi^  =  -\  for  3i<0  as  in  (1.6),  and  m^^O  otherwise.  For 
M  to  be  continuous  throughout  n,  it  remains  only  to  ensure  that  ^  is  continuous  at  t  =  «>. 
At  /  =  =o,  the  integrand  of  (1.5)  has  a  singularity  of  type  f^"  with  -y^  =  deg/7  + 
24»i  ('"i-Pi)  -  deg/7  +  (34_  +  2i*i,'»i-  Fo''  continuity  we  need  -yx<-l,  i-c,  degy7< 
-(l  +  3ij-  Ei.43-  Kife,  is  of  type  (1.6),  then  -(l  +  p^J  6  (-1,0]  and  -Si^i^m^-  ^-1, 
and  this  condition  becomes  dtgpsK-2.  On  the  other  hand  if  Jt^  is  not  of  type  (1.6),  then 
-(l  +  Pi )  e(-2, -1]  and  -2i?«i '"t"*^'  ^^^  again  we  have  degp^K-2.  This  completes 
the  proof  of  Theorem  3  except  for  the  assertion  regarding  maxima  and  minima. 

We  treat  the  case  of  a  maximum,  that  of  a  minimum  being  analogous.  Let  u  €  C(iT) 
be  a  nonconstant  solution  to  Problem  O.  By  the  maximum  principle  for  harmonic  func- 
tions, u  must  attain  its  maximum  m„^  somewhere  on  dCl.   Suppose  first  that  this  occurs  at 

a  point  C  ^  Ti  with  Ot  =  -j-,  or  at  a  vertex  adjacent  to  such  a  segment  r^ .  Then  u  ^  u^^ 
on  Ft  and  on  any  successive  adjacent  boundary  segments  F^  ^.j,  ^^ +2>  •••  ^^  ^ic~'T-  ^ 
the  entire  boundary  had  8i  =  -?-,  u  would  be  constant,  which  we  have  ruled  out.  There- 
fore there  must  be  a  final  boundary  segment  Ft +Ai  of  this  type,  and  the  vertex  z^+^+i  is 
then  a  point  with  w  =  «n^  at  which  (1.6)  holds,  as  asserted. 

Next,  suppose  u-u^^^  at  some  point  ?  ^  Fl  with  6i  <-j-.   By  the  representation  (1.5) 

just  established,  /  maps  a  segment  of  F^  around  I  onto  a  line  segment  oriented  at  angle 

0i  to  the  real  axis.    For  u  to  have  a  maximum  at  ?  we  must  have/' (0  =  0,  i.e.  p(t,)~0. 

But  in  this  case  (1.5)  implies  that  values  «>"„jaj,  will  be  attained  at  points  near  4  interior 
to  n.   This  is  a  contradiction. 

Finally,  suppose  u  attains  its  maximum  at  a  vertex  z^  with  e^-i.e^  <  -j-.  Again  we 
know  that  the  behavior  near  r^  is  given  by  (1.5),  so  that  /  maps  an  arc  of  dD,  around  z^ 
onto  two  straight  line  segments  with  a  comer  between  them.  If  «  is  to  have  a  maximum 
at  Z;,  /(z)  must  move  rightward  at  an  angle  8i_i  to  the  real  axis  as  z  moves  along  Vj^^^^  to 
Z;,  then  leftward  at  an  angle  6^-1-11  to  the  real  axis  as  z  continues  past  z^,  along  F^,,  and  the 
image  region  bounded  by  these  two  line  segments  must  he  in  the  half  plane  ReH'<Re/(z;). 
This  implies  0i<6i_,,  i.e.  (1.6),  as  claimed. 
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